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Modeling  of  the  Geosynchronous 
Orbit  Plasma  Environment -Part  I 


I.  INTRODUCTION 

Although  the  role  of  the  environment  In  generating  spacecraft  potential  var- 

i .4* 

iatlona  at  geosynchronous  orbit  is  well  documented,  variations  in  the  ambient 

environment  itself  have  not  been  well-defined.  Several  studies  have  been  made  of 

5-8 

the  geosynchronous  environment,  but  none  have  attempted  an  analytic  formu- 
lation of  Die  various  parameters  needed  by  Die  user  community  in  modeling  the 
spacecraft  charging  phenomenon.  This  initial  paper  describes  the  baste  set  of 
parameters  required  to  formulate  such  a model  of  Die  ambient  plasma  environ- 
ment and  outlines  a systematic  procedure  for  constructing  a simple  analytic  model 
that  includes  Uie  effect#  of  local  time  and  geomagnetic  activity.  Observational 
data  from  the  AT8-5  satellite  are  analysed  using  this  procedure  to  give  a prelim- 
inary analytic  description  of  Ote  geosynchronous  environment  in  the  form  of  a 
1-X)H  Til  A N program.  Although  not  intended  as  a detailed  description  of  Die  envi- 
ronment. Die  model  is  used  in  evaluating  local  time  variations  of  the  plasma  fol- 
lowing a plasma  injection  event  at  geosynchronous  orbit  (Uie  condition  most  likely 
to  foster  charging). 


CHecelved  for  publication  13  December  1977) 

#Due  to  the  number  of  references  to  be  included  as  footnotes  oa  this  page,  the 
reader  Is  referred  to  the  list  of  references,  page  35. 
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2.  PLASMA  CHARACTERIZATION 


The  first  step  in  constructing  any  model  is  the  determination  of  the  critical 
parameters  necessary  for  description  of  the  phenomenon.  In  plasma  physics, 
classical  studies  have  shown  that  to  describe  a plasma  completely,  the  electric 
field,  magnetic  field,  and  particle  distribution  functions  must  be  known.  Several 
models  of  the  magnetospheric  electric  and  magnetic  fields  exist  in  the  open 
literature  so  they  will  not  be  discussed  here.  This  paper  will  concentrate  instead 
on  the  so-called  plasma  distribution  function  - its  definition,  its  use,  and  its 
measurement. 

What  is  the  plasma  distribution  function?  Although  in  general  a complex 
concept,  the  distribution  function  is,  stripped  of  its  physical  and  mathematical 
niceties,  a function  which  describes  how  many  particles  exist  within  a tiny  volume 
of  space  with  a velocity  in  a certain  direction.  Turning  briefly  to  the  mathemat- 
ical details,  one  notes  that  F,  the  distribution  function,  is  given  by: 

F(X.  Y.  Z.  Vx.  /y.  Vz)  (1) 

such  that 

F(X,  Y,  Z.  Vx,  Vy,  V^)  dx  dy  da  dv^  dvy  dv^ 
is  the  number  of  particles  In  the  velocity  space 

dv^  dVy  dv^  and  spatial  volume  dx  dy  ds 


where 


X,  Y,  Z « spatial  coordinates 


and 


V . V8  ■ velocity  components. 

If  F is  integrated  over  oil  velocities  and  positions  within  a given  volume,  the  total 
number  of  particles  in  that  volume  is  obtained  (the  analysis  presented  in  this 
paper  assumes  a colltsionleM  plasma  with  rvo  particle  sources  or  sinks). 

For  the  purposes  of  this  study,  the  plasma  will  be  considered  to  be  Isotropic; 
that  is.  there  are  just  as  many  particles  traveling  in  direction  X as  In  direction  Y 
or  Z.  This  simplifies  the  analysis  greatly.  (Only  future  study  will  reveal  just 
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how  good  or  bad  such  an  assumption  is  and  this  factor  should  be  kept  in  mind  in 
the  following. ) Changing  to  spherical  coordinates,  one  notes: 

FiX,  Y,  Z,  Vx.  V V^)  dx  dy  dz  dv^  dv,_  dv 


y z ' x ~y  — z 
2 


a f(X.  Y,  Z,  v;  dx  dy  dz  sin  ev  d^  d<$v  dv 
a Y,  Z,  v)  dx  dy  dz  (4  rv2  dv) 


(2) 


where 


2 2 2 1/ ^ 


®v’  \ * anSular  coordinates  of  the  velocity  vector 
f « isotropic  distribution  function 


And  we  have  integrated  over  0 v and 


2*  * 


sin  fl  dO  « 4r 
v v 


(3) 


A commonly  encountered  distribution  function  is  (ho  so-called  Maxwellian 
distribution: 


/ m.  \3^2  -m.v2 /2kY 

«V«i(^)  • iJ  1 


(4) 


where 

• number  density  of  species  i 
»«t  • mass  of  species  i 
Tj  » temperature  of  species  I 
Vj  = velocity  of  l species 
k * Boltzmann  constant 
f * distribution  function  in  #ec3/baQ 
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Although  most  plasma  distributions  in  space  are  neither  Maxwellian  nor  isotropic, 
these  assumptions  are  commonly  made  in  characterizing  a plasma,  in  order  to 
reduce  the  number  of  parameters  necessary  for  description.  Further,  Eq.  (4) 
can  be  used  in  the  calculation  of  the  first  four  plasma  moments  that  will  generate 
the  model.  For  a Maxwellian  particle  distribution,  they  are: 


where 


s 

(Pj)  » number  density  for  species  l (number/ cm  ) 

(NFj)  **  number  flux  for  species  l (number/ cm*  sec-ar) 

, , . 2 .• 

(Pj)  • pressure  for  species  l (oynea/ cm  ) 

(EFj)  • energy  flux  for  species  i (ergs/ cm  sec-ar) 

Tim  use  of  moments  is  slmil  to  expanding  a function  in  a Taylor  series. 
Tim  momenta  are,  in  statistical  terms,  the  expectation  values  of  a variable,  that 
is,  the  average  value,  the  standard  deviation,  and  so  on).  For  a plasma,  the 
moments  of  the  velocity  arc  taken  as  follows:  (v >,  (v  ),  (v4) , and  (v  ).  In 
Eqa.  (5),  ($).  (?)  and  (8).  these  moments  have  been  multiplied  by  constants  to 
give  physically  meaningful  quantities.  For  example,  the  <v‘)  moment  has  been 
multiplied  by  4 r and  1/5  n»j  to  give  Uic  isotropic  pressure,  4*“  <1/3  v2}. 


is  two -thirds  of  the  so-called  energy  dessliy  which  is  in  unite  of 
energy  density  is  ihs  total  energy  per  unit  volume  of  the  plasma. 


erga/csr . 
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The  moments,  easily  derived  from  the  UCSD  ATS-5  data  (the  instrument 
employed  in  this  study),  can  be  used  to  derive,  a Maxwellian  and  a "2 -Maxwellian" 
fit  to  the  distribution  function,  as  will  be  demonstrated. 

An  interesting  and  important  feature  of  the  four  moments  is  that  they  (unlike 
the  plasma  temperature)  can  be  averaged  together  to  give  a phy  ically  meaningful 
average  quantity.  To  prove  this,  assume  a particle  distribution  is  given  by  f^. 
Then  the  moments  for  this  function  are  given  by 


M f ^ v2  dv 


where  M is  the  moment  to  be  found  for  distribution  i.  Now  assume  we  want  to 
find  (M)  for  a distribution  function  f where,  f is  the  average  of  several  distribu- 
tion functions  f,: 


N 


1*1 


Then  (M)  cat*  be  defined  as: 


<M>  • J 


Mfv  dv 


This  gives: 


(M)./.Mfv2dv«  f jM 

5 6 ' i«l  ' 

* k%(]  Wftv‘2  • iZ 

t«t\o  ! i"l 


t . 
v dv 


(MV 


Thus  groups  of  the  four  moments  can  each  be  averaged  to  give  a new  average  set 
of  moments  in  a i entirely  consistent  fashion.  These,  in  turn,  can  be  used  to 
derive  a new.  averaged  distribution  function. 


ii 


T-c-Si  «*»» 


tJ ^ ^.vi^v^igK 


r3.V'AJ  VJ^rWJ^VW’^1!**  ift'SSVS*  5. . 


In  Appendix  A,  it  is  demonstrated  that  if,  instead  of  Eq.  (1),  the  distribu- 
tion function  is  assumed  to  be  the  sum  of  2 Maxwellians  or  2 plasma  components 
for  a single  species  i,  then  n^,  n2i,  T^.  and  can  be  found  from  Eqs.  (5), 
(6),  (7),  and  (8)  such  that: 


-m,V2/2kT 
e 1 11  + 


2i 


-m.V' 


•/2kT, 


2i 


(9) 


In  Figure  1 are  plotted  actual  distribution  functions,  the  Maxwellian  fit,  and  the 
?■  Maxwellian  fit.  Neither  exactly  fits  the  distribution  functions  but,  of  the  two 
approximations,  the  2 Maxwellian  fit  gives  a much  better  representation  of  the 
actual  distribution  function  (for  most  studies  of  the  effects  of  the  ambient  environ- 
ment, a single  Maxwellian  has, been  considered  adequate  — clearly  a dubious 
assumption).  In  this  paper,  the  parameters  necessary  to  compute  both  the  single 
Maxwellian  and  2 Maxwellian  distributions  will  be  derived.  Figure  1,  however, 
should  be  remembered  in  considering  the  accuracy  of  either  representation. 

The  second  moment,  (NF^) , may  be  utilized  in  the  calculation  of  another 
important  quantity,  the  current  per  unit  area  to  the  spacecraft.  As  the  charge 
striking  a unit  area  of  the  spacecraft  per  unit  time  is  of  concern,  the  vector 
velocity  v.  relative  to  the  unit  normal  n to  the  surface  must  be  taken  into  account. 
Also,  only  particles  entering  one  half  of  the  sphere  (that  is,  particles  reaching 
the  satellite  surface  from  one  side  only)  are  considered.  Thus: 

V’l  /v5'1"3'’ 

ff/2  2ff  oo 

= q i J sin  9y  cos  0y  d9y  J d<fry  J f vJ  dv  (10) 

o oo 


\!/2 

2 \ rrm^  J 


<H  7 <NFt) 


where 

o-  = charge  on  species  (coulombs) 

J,  •»  current  por  unit  area  (amps/ cm  ). 
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Another  quantity  of  general  usage  is  the  mean  energy.  The  mean  energy, 

(Ej)  (that  is,  the  average  energy  per  particle),  for  a Maxwellian  distribution  is 
defined  as  the  ratio  of  the  energy  density  to  the  number  density  (where  (P^)  is  two- 
thirds  of  the  energy  density): 

(Et>  - | (P* >/<ni)  - f kT,  (11) 

where  kT.,  the  quantity  used  in  defining  the  Maxwellian  distribution,  can  be  ob- 
tained from  (E.)  by  means  of  thi3  relation. 

The  distribution  function  (approximated  either  by  a Maxwellian  or.  better, 

2 Maxwellian  fit),  the  mean  energy,  and  current  are  the  three  quantities  most 
often  required  by  programs  designed  to  calculate  spacecraft -plasma  interactions. 
Thus,  given  the  four  moments  and  Eqs.  (4),  (9),  (10),  and  (11),  a description  of 
the  plasma  environment  can  be  derived  which  meets  the  majority  of  the  user  com- 
munity's needs. 

3.  METHOD  OF  ANALYSIS 

In  the  preceding  section,  the  parameters  necessary  to  define  the  ambient 
plasma  were  described.  All  exhibit  large  temporal  and  spatial  variations  at  geo- 
synchronous orbit,  making  their  accurate  determination  quite  difficult.  In  fact, 
it  is  this  extreme  variation,  particularly  during  geomagnetic  storms,  that  has 
prevented  analytic  models  from  being  generated.  Tiro  method  to  be  outlined  is, 
as  a result,  a first-ordor  approximation.  Its  usefulness  can  be  determined  only 
in  retrospect  by  how  well  it  predicts  conditions  at  geosynchronous  orbit.  We  con, 
howover,  describe  *i>e  characteristics  of  a meaningful  model  and  build  those 
explicitly  Into  the  umilysls  process. 

What  characteristics  must  a model  of  the  geosynchronous  environment 
possess  ? In  particular-,  how  can  tiro  observed  temporal  and  spatial  variations  be 
accounted  for?  The  scale  of  the  variations  give  important  clues  to  tiro  answer. 
First,  changes  in  geomagnetic  activity  correlate  with  the  largest  temporal  varia- 
tions in  the  plasma.  Secondly,  spatial  variations  at  geosynchronous  orbit  trans- 
late Into  local  time  (or  position  relative  to  the  sun)  variations.  Thus,  the  two 
important  variations  are  geomagnetic  activity  and  local  time. 

Geomagnetic  activity  has  been  defined  historically  in  terms  of  geomagnetic 
storms  at  the  earth's  surface.  9 These  storms  appear  as  Initially  small  in- 
creases in  the  earth's  horizontal  magnetic  field  amplitude  <~10  - 100  y, 

9.  Rostoker.  G.  (1972)  Geomagnetic  indices.  Rev.  Gcophya.  j£(No.4):035. 
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where  a y is  10‘  of  a gauss  - the  earth's  field  is  -0. 3 gauss)  followed  by  a 
rapid  ~ 1000  y decrease.  This  decrease  may  last  a day  or  more  after  which  the 
field  slowly  (-week)  recovers  to  its  quiet  value.  Such  events  are  believed  to  be 
the  result  of  a compression  of  the  earth's  magnetic  field  by  the  solar  wind  and 
are  generally  accompanied  by  auroral  activity,  ionospheric  perturbations,  par- 
ticle fluxes  at  geosynchronous  orbit,  and  so  on.  A typical  magnetic  record  is 
shown  in  Figure  2. 
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Figure  2.  Kxample  of  a Geomagnetic  Storm  Showing  the  SSC  (storm  sudden 
commencement)  ari^  the  Resulting  Magnetic  Variations  at  the  Fur  til's  Surface 
(Akasofu  and  Cnapman.  1Q72) 


Traditionally,  the  level  of  > * -id-wide  geomag-et,c  activity  as  represented 
by  fluctuations  in  the  earth's  magnetic  field  has  been  defined  by  the  3 hour  K 
Index  which  goes  from  0 to  9 in  etc  ,.e  of  1/3.  This  is  a quasi -logartttimic  qur  ntlty 
whicn  is  supposed  tc  r^pros-nt  the  pooittve  and  negative  deviation  of  an  "average" 
magnetometer  at  mid-latitudes.  The  u^  index  la  directly  calculated  from  the  K 
index  and,  when  multipUed  by  2 y's  is  supposed  to  give  the  magnetic  deviation  In 
Y s.  Those  two  quantities,  duo  to  the  localized  nature  of  most  geomagnetic 
acUvlty,  arc  underestimates  of  ncaial  gecmegnetic  activity.  They  are,  however. 
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the  only  quantities  currently  available.  In  this  study,  A the  daily  sum  of  a , 

a r r 

will  be  employed  as  the  index  of  geomagnetic  activity. 

The  determination  of  how  the  distribution  function  changes  in  local  time  and 
with  Ap  is  a significantly  different  problem  from  determining  average  conditions 
or  extreme  values.  This  is  particularly  true  in  a study  of  charge  buildup  for  the 
proper  relationship  between  A^  and  local  time  relative  to  the  sun  of  the  satellite 
and  the  interrelationship  between  parameters  must  be  preserved  simultaneously. 
Finally,  in  order  to  make  the  model  of  general  utility,  it  should  be  expressible  in 
an  analytical  form.  To  fulfill  these  requirements,  a multiple  linear  regression 
analysis  was  carried  out. 

In  multiple  regression,  the  coefficients  of  a simple  analytic  expression  are 
adjusted  so  that  the  deviations  between  the  observed  and  predicted  values  as 
described  by  the  chi-squared  statistic  are  minimized.  For  this  initial  study,  a 
linear  relation  in  daily  and  a diurnal  pnd  semidiurnal  variation  in  local  time  t 
were  selected.  Thai:  is,  a collection  of  observations  varying  in  local  time  and 


of  a given 


moment  were  fit  using  the  following  equation: 


M(t.Ap)  = (aQ  + ajAp)(bo  + bj  cos  (t  + tj)  + b2  cos  (~r-(t  + tj))j  (12) 


where:  M(t,  Ap)  « predicted  value  of  the  moment  M at  local  time  t and  for  the 
daily  geomagnetic  activity  index  Ap,  and  aQ,  a^,  bQ,  b^.  b2,  t,  • coefficients 
determined  by  the  regression.  Using  a multiple  regression  progrom,  from 
Bevington. 10  we  note  that  aQ.  bQ,  bj.  b2,  tJ(  and  ^ were  determined  for  each 
set  of  moments,  A^,  and  t.  The  standard  deviation  for  each  point  was  assumed 
to  be  the  standard  deviation  of  each  set  of  momenta  for  analysis  purposes. 

In  summary,  the  plasma  environment  ut  geosynchronous  orbit  varies  with 
geomagnetic  activity  and  local  time  (that  is,  the  position  of  the  satellite  relative 
to  the  sun).  To  preserve  tills  relationship,  each  parameter  is  fit,  using  multiple 
linear  regression,  to  a simple  function  in  and  local  time.  Each  parameter, 
predicted  for  the  same  value  of  A^  and  local  time,  can  then  be  combined  to  give 
the  distribution  function  and  other  parameters  of  interest.  A FORTRAN  listing  of 
the  resulting  model  is  given  in  Appendix  A. 


10.  Bevington.  P.  R.  (196! 

Sciences,  McGrow - 

An  was  employed  In  Uou  of  su,  as  the  spread  In  sul In  the  data  base  was  inade- 
quate in  local  time  for  meaningful  comparisons.  The  use  of  A„  also  allows  Uto 
analytic  expression  to  include  the  observation  that  particle  Injections  at  geo- 
synchronous orbit  produce  effects  lasting  up  to  24  hours. 
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Data  Reduction  and  Error  Analysis  for  the  Physical 


4.  DATA  BASE 


To  test  the  method  just  described,  measurements  made  by  the  University  of 
California  at  San  Diego  (UCSD)  plasma  experiment  on  ATS-5  were  analyzed. 

In  August  1969,  ATS-5  was  launched  into  geosynchronous  orbit  (6.  6 RE>  near 
105°W.  The  satellite  spin  period  is  0.  79  sec  and  the  spin  axis  is  aligned  with 
that  of  the  earth's.  The  UCSD  instrument  consists  of  four  cylindrical  plate  spec- 
trometers: two  pairs  of  electron  and  positive  ion  detectors  directed  parallel  and 
perpendicular  to  the  spin  axis.  Three  simultaneous  measurements,  0.26  sec 
long,  can  be  taken  every  0. 32  sec.  A complete  spectrum  of  64  steps  of  energy 
channels  (two  are  background,  the  others  are  each  112  percent  of  the  previous 
channel  starting  at  51. 6 eV  and  ending  at  51.  6 keV)  can  be  taken  in  20.48  sec. 

Other  modes  are  possible,  but  for  this  study  only  two  to  four  minute  averages  of 
the  four  moments  of  the  distribution  function  are  studied.  A more  detailed  descrip- 
tion of  the  calculation  of  the  four  moments  from  the  UCSD  data  is  given  in  Appen  * 
a ix  A. 

As  the  primary  purpose  of  tills  paper  is  to  outline  a procedure  for  obtaining 
an  analytic  formulation  of  die  geosynchronous  plasma  useful  in  studing  charging 
effects,  we  do  not  use  a large  data  set  nor  are  the  days  selected  at  random  (see 
the  following).  Ten  days  of  hourly  measurements  wore  chosen  which  covered  a 
wide  range  of  geomagnetic  activity.  Data  gaps  and  other  singularities  were  Ig- 
nored and  interpolated  values  employed  (even  so.  this  involved  only  about  0 of  the 
240  observations).  Table  1 lists  the  ton  days  chosen  and  the  corresponding  3 -hour 
ap  and  dally  Ap  values  (see  Rostokor,  9 for  a review  of  geomagnetic  indices). 
Initially  250  sets  of  moments  were  selected  (one  set  every  hour  on  the  half-hour 
plus  the  values  at  plasma  injection,  if  clearly  identifiable),  17»e  plasma  compo- 
nents perpendicular  and  parallel  to  the  satellite  spin  axis  were  averaged  together, 
assuming  equal  weights.  Finally,  to  facilitate  a comparison  with  the  3-hour  s^ 
index,  3 -hour  averages  were  taken.  It  was  tills  data  base  of  80  values  per  moment 
which  was  analyzed. 

Hie  data  base  was  not  selected  at  random  to  avoid  errors  In  the  local  time 

variations  which  are  affected  by  single  injection  events.  As  discussed  In 

5 4 

Deforest  and  Mcllwuln  or  Garrett  et  ai,  on  lt\jection  event  is  the  sudden  appear- 
ance of  hot,  plasma  sheet  plasma  near  local  midnight  at  geosynchronous  orbit. 
Unfortunately,  these  events  are  by  their  very  nature  random  in  occurrence.  Such 
events  cause  order  of  magnitude  changes  in  the  plasma  conditions  at  geosynchron- 
ous orbit  and  tend  to  skew  any  gross  averages  of  the  geosynchronous  environment, 

*The  electron  data  for  1972,  by  which  time  these  detectors  had  degraded,  were 
corrected  using  a factor  of  48  for  the  parallel  detector  and  0 for  foe  detector 
perpendicular  to  foe  spin  axis. 


Table  1.  Geomagnetic  Activity  for  Days  Analyzed 


Year 

Day 

Hourly  AP 

Daily  AP 

1972 

217 

67 

236 

132 

27 

56 

27 

111 

400 

1056 

1970 

348 

22 

80 

236 

48 

32 

15 

32 

56 

521 

1972 

223 

9 

6 

12 

18 

32 

27 

15 

27 

146 

1969 

326 

15 

18 

12 

7 

5 

5 

12 

6 

80 

1970 

273 

6 

18 

18 

5 

7 

6 

4 

9 

73 

1970 

272 

3 

5 

6 

5 

2 

3 

2 

12 

38 

1969 

209 

3 

4 

9 

3 

4 

2 

5 

3 

33 

1970 

25 

4 

3 

2 

0 

2 

2 

3 

4 

20 

1971 

87 

5 

2 

0 

0 

2 

2 

4 

3 

18 

1970 

345 

0 

2 

3 

3 

2 

0 

2 

2 
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As  the  purpose  of  this  analysis  is  to  model  the  time-history  of  a charging  event, 
we  have  selected  days  exhibiting  single  injection  events  near  satellite  midnight, 
for  the  lower  levels  of  geomagnetic  activity.  Such  steps  were  based  on  the  obser- 
vation that  spacecraft  charging  has  been  found  to  occur  primarily  in  conjunction 
with  Individual  plasma  injections  near  midnight.  This  bias  should  be  considered 
when  evaluating  the  generality  of  the  modal. 


S.  RESULTS 

A variety  of  parameters  and  combinations  of  parameters  were  calculated 
for  this  study.  Considering  the  sparsencss  of  the  data  base,  we  have  limited  the 
study  to  a formulation  of  the  four  momenta  in  terms  of  linear  regression -derived 
analytic  expressions.  The  currents  and  the  parameters  necesssry  for  a 
Maxwellian  and  a 2 Maxwellian  fit  to  the  distribution  function  are  computed  horn 
these  predicted  quantities.  As  a result  of  the  simplicity  of  the  fitted  function  (in 
particular  the  assumed  linearity  of  A^)  and  the  paucity  of  data,  some  negative 
values  result  for  the  predicted  values.  As  these,  in  turn,  result  in  fictitious 
values  for  the  derived  quantities,  in  the  FORTRAN  formulation  of  the  problem, 
these  values  have  been  corrected  by  estimates  (sec  Table  2).  With  this  caveat 

*The  average  energy  was  calculated  from  the  data  and  was  derived  from  the  fitted 
values.  The  latter  values  are  used  in  tilts  preliminary  model  to  pi'eserve  the 
internal  consistency  of  the  model,  since  the  difference#  between  the  values  so 
derived  are  within  the  accuracy  of  the  model. 


Table  2.  Default  Values  and  Standard  Deviations  for  Model 


" — m 

V ariable 

Units 

Electrons 

Std.  Dev.  /Default 

Ions 

Std.  Dev.  /Default 

Density 

no#  / cm  9 

0.  67 

0.02 

0.  55 

0.33 

Pressure 

o 

dynes/ cm 

2.6  X 10*9 

4 X 10'11 

4.  9 X 10~9 

4 X 19~9 

Energy  flux 

erg/(cm2-sec 

-sr) 

2.  1 

0.  08 

0.  13 

0.  09 

Number  flux 

2 

no#/ (cm  sec 

-sr) 

1.4  X 108 

4 X 106 

4 X 106 

3 X 106 

Mean  Energy 

eV 

(1600) 

1000 

(1820) 

Current 

nAraps/cm2 

0.  07 

0.  002 

0.  002 

0.  0015 

N1 

no#/cm9 

2 

N2 

no#/cm2 

0.  04 

T1 

eV 

250 

T2 

eV 

20000 

in  mind,  the  model  predictions  will  be  discussed  and,  where  possible,  compared 
with  other  observations. 

Tables  3 through  22  list,  by  parameter,  the  values  predicted  by  the  model 
(sec  FORTRAN  listing  in  Appendices).  The  parameters  are  given  as  a (unction 
of  the  A,  vulues  corresponding  to  overage  daily  values  of  K of  0 , 1 , 2 , 3 , 4 , 

1}  p U v U U M 

5 , 6„,  7,8,  and  8 . As  the  data  set  of  A , values  Includes  only  values  as  high 
o*  o o o o p ^ 

as  1056  (that  is,  K ■ 7^),  the  values  for  1656  and  3200  should  be  treated  as 
extrapolations.  For  each  of  the  four  moments  and  tor  the  current,  die  standard 
deviation  associated  wttli  the  fit  to  the  actual  data  is  given.  This  is  indicative  of 
the  error  that  should  bo  assigned  each  predicted  value.  As  would  be  expected  for 
the  small  data  set  employed  in  tills  study,  the  predicted  errors  are  a large  fraction 
of  the  model  values.  Also,  the  environment  actually  exhibits  large  variations  much 
as  these. 

In  consideration  of  the  errors  associated  with  the  four  moments,  it  is  not  jus- 
tified to  draw  more  than  a few  generalized  conclusions  from  the  predictions  of  ihto 
preliminary  model.  In  general,  the  ranges  of  die  four  moments  wall  approximate 
those  given  by  DeForest  and  Mdlwatn5  in  their  Table  i if  their  "maximum1'  values 
arc  taken  as  corresponding  to  an  A^  of  105(5  (average  1C^  of  ?fi),  their  "typical" 
values  as  corresponding  to  Ap  of  - 120  (average  of  3fl},  end  "tiilnlmuto1"  to  an 
of  -0  (average  of  0Q).  The  agreement  is  excellent  for  the  ions  os  to  actual 
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Table  6-  Ion  Pressure  X 1010  (dynes/  cm^) 
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amplitude  and  locel  time.  The  agreement  is  also  quite  good  for  the  energy  flux 
and  number  flux  'or  tne  electrons,  but  not  as  good  with  the  density  and  pressure 
(the  model  predicts  two  peaks  - one  near  midnight  and  the  other  near  0930  in  con- 
trast to  their  single  peak  near  midnight).  In  consideration  of  the  greater  varia- 
bility of  the  electrons,  this  discrepancy  is  probably  a real  deviation. 

Although  skewed  somewhat  by  the  problems  discussed  in  earlier  sections  of 
the  paper  (that  is,  injections  and  a less  then  random  selection  of  events),  the 
amplitudes  of  the  four  moments  predicted  by  the  model  for  the  electrons  and  ions 
demonstrate  consistent  local  time  variations.  For  example,  all  show  a sharp 
minimum  between  1630  to  1930  local  time  for  moderate  to  high  le\els  of  geomag- 
netic activity  - the  minima  being  on  the  order  of  50  percent  of  the  maximum 
values.  At  low  values  of  geomagnetic  activity  and  for  the  derived  quantities,  thiB 
trend  is  lost  or  greatly  reduced.  The  minimum  is  obvious  in  most  spectrograms 
returned  by  ATS-5  and  ATS-6  and  probably  reflects  the  sharp  edge  of  the  plasma 
sheet/ injection  boundary  in  the  evening  hours  reported  by  many  others1*'  ‘ 
which  marks  the  boundary  between  eastward  drifting  low  energy  ions  and  westward 
drifting  high  energy  ions. 

In  contrast  to  the  minima,  the  peaks  at  high  geomagnetic  activity  appear 
broad  in  extent  and  well-defined.  At  low  levels  of  A , all  moments  peak  between 
0130  and  0430.  The  mean  energies  for  the  electrons  are  about  one -fourth  of 
expected  valuta  and  show  no  strong  local  time  or  A^  variation,  although  there  may 
be  a slight  tendency  for  higher  energies  to  occur,  on  the  average  near  0430  local 
time.  The  mean  ion  energies,  in  comparison,  peak  near  midnight  moving  to  ear- 
lier hours  (~2000)  as  the  geomagnetic  activity  decreases.  Again,  tills  would  be 
consistent  with  the  movement  of  the  hot  plasma  sheet/ injection  boundary  toward 
die  evening  hours  at  geosynchronous  orbit  during  geomagnetic  activity. 


The  derived  quantities  (mean  energy,  Nj. 
tamlnutod  by  minimum  value  estimates  below  A 


N 


3>  , and  Tg),  ulthough  con- 

- 56,  reveal  several  features. 


As  j . it  described,  the  election  mean  energy  shows  no  strong  variation  with  local 
time  or  A . Hi  is  is  consistent  with  the  observations  of  ATS-5  data  reported  by 


11.  Vasyliur.as,  V.M.  (1968)  A survey  of  low  energy  elections  In  the  evening 
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13.  Mauk.  B. , and  Mcllwain,  C.E.  (1974)  Correlation  of  K_  with  the  substorm - 

Injected  plasma  boundary,  J.  Geopltya.  Res.  79:3103. 
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14 

Inouye.  The  estimates  of  the  standard  deviations  of  these  values  are  on  the 
order  of  ~ 1500  eV  for  the  electrons  and  ~2000  eV  for  the  ions,  quantities  of  sim- 
ilar magnitude  for  the  electrons  as  the  change  in  the  predicted  values.  This 
trend  is  evident  in.  the  2 Maxwellian  fits  exci.pt,  interestingly,  the  high  energy 
temperature  component  of  the  electrons  and  ions.  For  them  a trend  toward 
higher  values  for  the  daily  average  temperature  at  lower  levels  of  geomagnetic 
activity  is  significant  and  and  in  contrast  to  the  increase  in  the  daily  average  of 
all  other  parameters.  This  would  lend  support  to  the  assumption  of  Inouye14  and 
Stevens  et  al  that  the  plasma  temperature  goes  up  as  the  current  goes  down. 

The  results  imply,  however,  that  the  effect  is  dependent  on  how  one  defines  the 
energy  or  temperature.  The  mean  energy  as  defined  in  this  study  is  — l/4  of  its 
expected  value  and  usually  increases.  In  fact,  the  temperature  as  represented  by 

T,  for  the  electrons  and  ions,  not  the  mean  energy,  better  represents  the  maxi- 

£ 7 4 15 

mum  values  estimated  by  other  studies.  ’ ’ Also,  the  maximum  temperature 

and  mean  energy  for  the  electrons  are  a factor  2 higher  if  maximum  instead  of 

average  values  are  used  in  the  study. 

On  the  assumption  that  the  2 Maxwellian  fit  is  a reasonable  approximation 
to  the  actual  plasma  distribution,  further  observations  can  be  made  based  on  the 
model:  First,  It  is  the  low  temperature  component  (Tj)  of  the  ions  and  electrons 
which  shows  the  largest  percentage  increase  as  geomagnetic  activity  increases. 
Related  to  this  is  the  observation  that  although  the  densities  Nj  and  Ng  increase 
markedly  with  geomagnetic  activity,  the  ratio,  Ng  / N j of  the  daily  averages  of  the 
densities  remain  roughly  equal  in  spite  of  largo  local  time  variations  in  the  ratios. 
This  may  Indicate  that  the  percentage  of  particles,  as  represented  by  their  ratio 
Ng/Nj,  in  the  two  populations  changes  little  on  a dally  basis,  whereas  their  tem- 
peratures Tj  and  T^  are  varying  Independent  of  each  other  (in  fact,  the  high  tem- 
perature components  may  be  decreasing  with  geomagnetic  activity).  In  considera- 
tion of  the  paucity  of  data  in  the  analysis,  however,  these  conclusions  are  tentative. 


b.  vtOUM.  tiSAGK 

In  this  section,  two  examples  of  model  usage  In  estimating  the  effocts  of 
charging  will  be  given.  The  model  wlU  be  employed  In  estimating  tho  environment 


14.  Inouye.  G.  Y,  (1976)  Spacecraft  potentials  in  a substorm  environment,  AIAA 
Progress  in  A stixjnaul'lcs'unQ  Aeronautics  Series,  Vol.  45,  pp.  103-120. 

’5.  DeForcst,  S.  E.  (1977)  Final  Report  for  i June  1976-30  November  1976, 
AFGL'TR  -77-0031.  
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conducive  to  charging  at  geosynchronous  orbit.  This  environment  can,  in  turn, 
be  joined  with  a spacecraft  charging  model  to  give  satellite  potentials. 

Of  general  interest  is  the  prediction  of  the  charging  environment  a satellite 
is  likely  to  encounter  during  a typical  mission.  The  basic  idea  is  to  predict  the 
expected  occurrence  frequency  of  during  the  mission.  This  has  been  done  for 
a ''typical"  time  period  in  Figure  3 where  a histogram  of  over  30  years  of  A 
values  has  been  plotted.  The  results  of  this  figure  can  be  utilized  in  two  ways. 
First,  the  probability  of  observing  a given  interval  of  A or,  equivalently,  a given 

r 

level  of  geomagnetic  activity  can  be  found;  for  example,  the  probability  of  observ- 
ing the  Ap  interval  56-120  is  33  percent  in  a given  time  period.  Similarly,  the 
probability  that  a value  of  A^  of  120  will  be  exceeded  during  a mission  is  given  by; 

P = 100%  - 14% (A  <32)  - 23%(32  < A :£  56) 

P P 

- 33% (5 6 < A £ 120)  ° 30%  . 

P 

The  value  of  A^  (120  in  this  example)  can  then  be  substituted  into  the  model 
(see  Appendices  for  a FORTRAN  listing  of  the  model)  along  with  a local  time  be- 
tween 0000  and  2400  to  give  the  ambient  environment  expected  33  percent  of  the 
time,  or  not  to  be  "exceeded"  more  than  30  percent  of  the  time. 

Substituting  the  A value  of  120  and  a LT  of  0130,  the  FORTRAN  subroutine 
* P 

returns,  one  notes: 


Number  density  (no#/cin  ) 

2 

Pressure  (dynes/ cm  ) 

2 

Energy  flux  (erg/cm  soc-sr) 
2 

Number  flux  (no#/cm  sec-sr) 

Mean  energy  3/2  kT  (eV) 

2 

Current  (amps /cm  ) 


Electron 

Population 

Ion 

Population 

1.49 

1.38 

3.31 X 10'9 

1,24  X 10'8 

1.09 

0.30 

2.57  x 10° 

107 

2100 

8400 

0.  13 

0.  005 

The  model  is  not  intended  to  be  used  to  give  total  average  plasma  "dosages" 
since  plasma  Injections  occurring  when  the  satellite  was  near  local  midnight 
were  preferentially  selectod.  Rather,  it  is  intended  to  give  the  environment 
most  conducive  to  charging  - that  is,  following  injections  when  the  satellite  is 
near  local  midnight. 

^ s 

Note;  ITiese  values  are  in  very  good  agreement  with  DeForost  and  Mcllwain  and 
Garrett  et  al4  for  "typical"  values  if  3/2  (T2)  is  the  "moan  energy.  " 


2 Maxwellian: 

Electron 

Population 

Ion 

Population 

N1  (no#/cm3) 

1.  13 

0.  54 

Tl  (eV) 

496 

21 

N2  (no# /cm3) 

0.  35 

0.84 

T2  (eV) 

4270 

9200 

The  approximate  distribution  function  (or  particle  spectrum)  can  then  bo  derived 
by  either  of  the  following: 

Electrons 


fQ(E)  «■  27.2  No 


(12a) 


30 


i 


Ions 


fjCE)  « 2 


B ( Ti  \'3/2  ‘E/TI 
. 14  X 10  Nj  ( 1000  / e 


(12b) 


where 


3 «*  ft 

fe,fj  = distribution  functions  (sec0  km*  ) 


N 3 1. 49  cm 
e 

Nj  ° 1.  38  cm 


-3 


-3 


= 2/3  (2100  eV)  = 1400  eV 
3 2/3  (8400  eV)  = 5600  eV 
* particle  energy  (eV) 


or: 


f'(E)  « 27 
e 


/ / 


Tl  \*3^2  -E/Tl 
1000/ 


e 


o 


+ N2 


T2  \*3/2  -E/T2 


(IM 

\ 1000/ 


(13a) 


-3/2 


f^(E)  - 2.  14  X 10°  \ Nl 


/Tl.  V°/fi  -E/Tl. 

i \Tooo/  e + N2 


T2t  \*3/2  -E/T2j. 


(Hi 

\100( 


i v looo  y e 


(13b) 


I 


where 


f‘ , fj  » distribution 

C 1 

Nl  « 1. 13  cm'3 
e 

Tlfi  » 496  eV 

N2„  » 0.35  cm'3 
e 

T2„  - 4270  eV 
e 


3 

functions  (sec0  km  ) 
N lj  ■ 0.  54  era  3 
Tlj  ■ 2 1 eV 
N2j  ■ 0.  84  cm  3 
T2j  ■ 9200  oV 


Tlie  resulting  distribution  functions  are  plotted  in  Figures  4 and  5.  It  is  recom- 
mended that  Eq.  (13)  be  used  for  Ap  greater  than  50,  and  Eq.  (12)  for  lower 
values. 


A second  use  of  the  model  is  in  making  real  time  estimates  of  the  plasma 
conditions.  A provisional  a^  value  and  a predicted  a^  value  for  the  next  3-hour 
period  are  currently  available  from  Air  Force  sources.  If  a running  sum  of  8 a^'s 
is  maintained,  the  model  can  be  used  to  predict  the  ambient  conditions  (that  is,  if 
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Figure  6.  Local  Time  Plot  of  the  Values  of  the  Four  Elec- 
tron Moments  Predicted  by  the  Model  for  A of  120 

P 


Figure  7.  Local  Time  Plot  of  the  Value#  of  the  Four  lou 
Moment#  Predicted  by  the  Model  for  Ap  of  12 0 
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the  preceding  7 ap  values  totaled  105  and  the  predicted  ap  was  15,  then  Ap  pre- 
dicted would  be  12  0 and  the  environment  versus  local  time  for  the  next  three 
hours  would  be  as  shown  in  Figures  6 and  7).  Likewise,  if  the  current  Ap  value 
is  120,  Figures  6 and  7 illustrate  the  average  conditions  following  an  injection 
near  the  local  midnight  portion  of  a satellite’s  orbit,  the  condition  most  likely  to 
foster  charging.  Extreme  conditions,  suitable  for  providing  a spacecraft  charging 
alert,  are  estimated  by  taking  the  largest  observed  or  predicted  a value  in  the 
desired  interval  and  multiplying  by  8 (if  the  highest  ap  was  15,  then  Ap  predicted 
would  be  120).  The  values  so  derived  can  then  be  inserted  in  a program  that 
computes  spacecraft  potential  to  give  an  "alert"  bulletin  versus  local  time. 


7.  CONCLUSION 

In  review,  we  have  outlined  a procedure  for  generating  an  analytic  formu- 
lation of  the  various  parameters  needed  by  researchers  seeking  to  model  the  geo- 
synchronous environment  and  the  interaction  of  that  environment  with  a space- 
craft. An  environmental  model  based  on  a limited  data  sot  (10  days)  was  analyzed 
under  this  procedure.  The  results  were  compared  with  other  observations  of  the 
geosynchronous  plasma.  Although  the  model  was  designed  to  analyze  plasma 
variations  following  a plasma  injection  near  the  midnight  portion  of  the  satellite 
orbit,  excluding  anisotropic  fluxes  and  orbital  effects,  it  included  geomagnetic 
and  local  time  variations.  Key  features  of  the  magnetosphere  such  as  the  plasma 
decrease  near  evening  were  reproduced,  and  various  trends  in  the  data  which  may 
be  significant  noted.  In  view  of  the  assumptions  and  size  of  the  data  base,  this 
model  is  considered  to  be  a preliminary  rather  than  a definitive  description  of  the 
ambient  geosynchronous  environment.  As  the  power  of  the  technique  has  been 
demonstrated,  it  Is  planned  to  extend  It  in  the  near  future  to  a much  more  compre- 
hensive data  base  and,  In  turn,  generate  a more  complete  model. 
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Appendix  A 

Calculation  of  Momoof* 


The  ATS-5  UCSD  instrument  makes  measurements  only  between  Si.  6 oV  and 
51.6  keV.  Further,  all  positive  tons  are  assumed  to  be  ionised  hydrogen  as  it  is 
not  possible  to  calculate  accurately  the  composition  of  the  ion  population.  Thus, 
the  four  moments  calculated  from  the  ATS-5  data  are  in  actuality  valid  only  for  the 
range  5 1. 6 eV  to  5 1.  6 keV  for  electrons  and  protons.  Do  Forest1  notes  that  the 
main  error  introduced  by  these  effects  Is  an  underestimation  of  the  particle  densi- 
ty and  a lack  of  equality  between  the  proton  and  electron  densities  as  a large  part 
of  the  particle  population  is  below  SO  eV.  With  this  in  mind,  the  formulas  for  the 
moments  in  Eqa.  (5),  (6).  (7),  and  (8)  are  approximated  for  the  electrons  by: 


(1A) 


(npv-  2 

E, 


d£ 


i x 


(2A) 
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(3A) 


(4A) 


where 


" summation  over  the  energies  from  51.  6 eV  to  51. 8 keV 


rn 

e 


* mass  of  electron 

3 enerffy  nux  (count  rate  divided  by  4.  3 X 10'5  cm2  sr 

for  ATS-5; 

« o.  12  for  ATS-5 


Hie  (on  moments  are  similarly  derived  by  replacing  m by  m..  the  mass  of 
a proton.  0 ' 
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Appendix  B 

Calculation  of  tf\e  2 Monwollion  DiUHbotioo 


The  distribution  (unction  ( is  of  central  Importance  in  various  schemes  for 
Uto  calculation  of  spacecraft  potential.  As  Figure  1 clearly  demonstrates,  a 
single  Maxwellian  distribution  is  often  inadequate  in  describing  Uto  actual  data. 
AlUtough  a 2 Maxwell  fit  is  also  only  an  approximation,  it  does  double  the  knowl- 
edge of  the  actual  distributions.  Further,  It  tends  to  divide  the  particle  distribu- 
tion Into  two  components:  a high  temperature,  relatively  low  density  component 
and  a low  temperature,  high  density  component  which  is  consistent  with  various 
observations  of  the  actual  stale  of  the  magnelospherlc  plasma.  For  these  reasons, 
to  Include  the  calculation  of  the  two  densities  (N ^ and  N.>)  and  temperatures  <Tj 
and  1%)  in  U»c  model  was  warranted. 

Assuming  a 2 Maxwellian  distribution  to  be  given  by  Eq.  (9)  and  making  Uto 
additional  assumption  that  the  four  moments  given  in  Appendix  A represent  ade- 
quate approximation  to  Ute  actual  distribution  function,  we  obtain: 


(SA> 

NiXl  + *WC2 

<6A) 

N1X1  + *VN2  * C3 

(7A) 

"lXI**V4-C4 

(8A) 
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where 


N ,N2  3 number  densities  for  species  i 

X1 

1/2 

= T £ , temperature  1 for  species  i 

X2 

1?2 

= Tg!  , temperature  2 for  species  i 

C1 

= <ni> 

S 

= (NFt) 

C3 

= (Pt)/K 

C4 

= <EF.) 

Solving: 

Niy 

C1‘N2 

(9A) 

N2  = 

C2'X1 

X2'X1 

(10A) 

X2 

C3-C2X1 

(11A) 

C2-C1X1 

(C2  ' 

■ CjCg)  x2  4 (CjC4  - c2c3)  Xj  + (C2  - c2c4)  = 0 

(12A) 

Equation  (12A)ts  a quadratic  equation  and  has  two  roots  of  the  form: 

V"  r- 

-B  + Vb2  - 4AC 

(13A) 

X1 

2A 

-B  - Vb2  - 4AC 
X1  „ — 


where 


A . C2  - C.  C, 


"1 


fc  r 


This  situation  may  appear  ambiguous,  but  it  is  obvious  by  the  symmetry  of 
Eqs.  (5A),  (6A),  (7A),  and  (8A)  that  if  we  choose  the  positive  sign  for  X^,  then 
X2  must  correspond  to  the  negative  sign.  The  actual  problem  is  that  it  is  possi- 
ble that  the  distribution  function  is  a single  Maxwellian  in  which  case  X„  (or  by 
symmetry,  Xj)  will  approach  infinity  and  the  number  density,  N„,  zero.  Some 
care  must  be  exercised  in  studying  trends  as  a result  of  this  effect. 

For  actual  data,  it  is  not  likely  that  imaginary  or  negative  values  will  be 
encountered.  For  values  derived  from  our  simple  model,  this  is  not  true:  The 
fitted  moments  may  be  negative  (in  which  case  imaginary  values  are  obtained),  or 
the  minimum  correction  values  may  result  in  negative  densities.  Fortunately, 
this  happened  for  only  7 of  the  80  values  calculated  for  the  electrons;  it  did  not 
occur  for  the  ions.  Default  values  representative  of  the  expected  minimum  values 
are  returned  by  the  program  in  these  cases  (see  Appendix  C).  Finally,  T^2  or 
T^2  may  be  negative  and  should  be  checked  in  any  general  application. 


Appendix  C 

FORTRAN  Listing  of  Environmental  Modal 

The  model  outlined  in  this  study  consists  of  two  subroutines:  MODEL  and 
MAXW.  These  two  programs  have  separate  purposes  and  can,  with  minor  changes, 
be  used  independently. 

Subroutine  MODEL  requires  as  input  the  daily  index  and  the  local  time 
LT  (a  real  number).  It  returns  a vector  XX  which  has  the  20  components: 

XX(1)  “ Electron  density  x 100  (number/cm3) 

XX(2)  ° Ion  density  x 100  (number/cm3) 

in  9 

XX(3)  ■ Electron  pressure  x 10  (dynes/cm6) 

XX(4)  ■ Ion  pressure  x 1010  (dynes/cm2) 

XX(5)  « Electron  energy  flux  x 100  (erg/ cm2  sec-sr) 

XX(6)  ■ Ion  energy  flux  x 100  (erg/ cm2  sec-sr) 

XX(7)  e Electron  number  flux  x 10’6  (number/cm2  sec-sr) 

XX(0)  * Ion  number  flux  x 10  (number/ cm  sec-sr) 

XX (0)  “ Electron  mean  energy  (eV) 

XX(io)  * Ion  mean  energy  (oV)* 

Moan  energy,  which  la  3/2  KT  for  a Maxwellian,  should  not  be  confused  with  the 
temperature,  KT,  used  in  generating  the  Maxwellian  or  2 Maxwellian  distribution. 
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4 2 

XX(ll)  = Electron  current  x 10  (n  amps/ cm  ) 

4 9. 

XX(12)  = Ion  current  x 10  (n  amps/cm  ') 

O 

XX(13)  = Electron  density  1 x 100  (number/cm  ) 

XX(14)  = Ion  density  1 x 100  (number/ cm"*) 

XX(15)  = Electron  temperature  1 (eV) 

XX(16)  = Ion  temperature  1 (eV) 

3 

XX(l7)  = Electron  density  2 x 100  (number/ cm  ) 

q 

XX(18)  = Ion  density  2 x 100  (number/cra  ) 

XX(19)  = Electron  temperature  2 (eV) 

XX(20)  » Ion  temperature  2 (eV) 

Subroutine  MODEL  requires  MAXW  but  it  can  easily  be  deleted  if  the 
2 Maxwellian  distribution  is  not  required.  Subroutine  MAXW  requires  a set  of 
four  moments: 

RHO  = density  (number/ cm  ) 

2 

FNO  » number  flux  (number/ cm  -sec-sr) 

n 

PR  = pressure  (dynes/ era  ) 

2 

FEN  » energy  flux  (erg/cm  -sec-sr) 

A value  of  1 is  used  for  electrons,  2 for  ions.  It  returns: 

O 

R1  ® density  1 (number/ cm  ) 

Y •>  temperature  1 (eV) 

R2  ■ density  2 (number/  cm  ) 

T « temperature  2 (eV) 

These  values  can  be  used  in  Ec.  (9)  to  give  an  approximation  to  the  distribution 
function. 

Table  2 of  the  main  report  gives  default  values  returned  by  the  program. 
These  are  the  most  reasonable  estimates  of  these  values  that  could  be  determined. 
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SU3P0UTIME  MOOEL<**>,LT,XX) 


REM.  Lf 

-OftT-*-  (( X+ir-W-yf*!^)  »F 

1 .3835c*02,-.4166E+02,  .2184F+02, 

Hr 


210it+01 

TAa>9F»'48T~r-7337e-84,-,49Sa£-0-lT-«80S8E-01 


, -» ‘’312F*  02, 

, , 4438P-01, 
,-.3287E+0i, 
r-rF«3£4ir 


1 . 9853c+02,-.27?0F+U2*  .2276E+01, -.2758E+C1 

-t — .2880 £+00,  r9757E~a-2,— .^543E-01tt~.4434E-01 
1 « 5195E+J 1,  -• 10  36E+02,  ,6420E*dl,  .2&09E+01 


, -.477*7*  01, 


1 »7684E+3 2, -• 4464E+01,  .1610E+0 2, 80 89E+00,-,3< 10r*01, 

■4 .464»S+«4t“.  4948e-»ai,-.'41»9&>»4-r^W7»E-wl7  -rg-369F--01-, 

1 .3579E*02,-.7375E+d2,  ,50‘»6E+Q2,  .2135E+Q2 ,-.10 7*E+02» 


1 .1858E+02,-. 5i58E+0a,  .3i09E+ai,-.6267E-Ul,-.*fl42E«-  00, 

-i — v 377&-E-0 1,-»  1237E- 01,-  .H+90E-81,-ui782E-O4,  .3-398€~02, 
i .4719E+0 2, -.6769E+02,  .3818E+Q2,  .7563F*ai,~.  M+2PF*02, 
■4 — rfrt4^+04,--=-a4-g5E  + fr8,-. »471E-fr*r-i.-ld»9E»98 , ,£WH^r 
1 .6232E+0 1,  -«  9452E*  00,  .i427E*Jl»-.i2d2£+0O,-.3"29F*OO, 

DATA  ( XX I m)  ,1=1,  8)  /2.U,  33  40.  ,8., 9. ,4., 3./ 

—Ft:  ASM 

00  2 1=1,8 

...  -XX  f4»-3 


00  1 J=+,9 
IPA-J/5 


ipt*j-5*t»a 


— ¥*4 

XI*  (LT 4-6,  5)  /3, 

— *PI-6. 28318 

IF(IPr.E?.l> Y=00S<Xi*TPI/8.) 

— 3F  ( It>Ta £9,  2 >■  V SIN (XI «TFiy8.0 

IP  (IBT.  EO.  3)  Y=  COS  tXl*  T*I/4.  ) 

■ — If^IPT,efi,4»Y~SIN<Xl«TPI7W)~  

J1=J+1 

4 XX<l)*(AP»*IF4)»»*)P(41,?)-»XXat 

XXtI)*XX(mX(l»I) 

TFOtKtl  >.LT,XMINtn)-FE»F=l  — 

2 IF (XX (I) . ir.XMIN(I) ) XX  <Y)  = XHIN(T) 

— XX  < 9>«XX(3)»9340./XX4H 

XX(10)=XX(4>*9360./XX(2» 

— XX  <44H-X*  PFP»l-«-6*3,  4415  9 

XX <12>*XX< 8) *1.6*3.14159 

— &e-a-£»4Ta 

K=I-1 

— RMO*XXt<n>  7100. 

FNO=XX(K*7>*1,E*06 

— pq»XXH-,4).-/+T»-fr+-4< 

FEN=XX<<*5» '100. 

— '‘■ALL  4AXX(R40,  FNOyWt,FeX,>>B,C,t),I> 
XX (K*13> ***100. 

— K-X'fK+lrS— 3 

X«(XM7M0*10g  , 

■3-XX1X+49V-3 

IFfXXttS*  ,t.r.80a,)G0  TO  5 
— IFtFi  49  »C3r3n,f-80-fO"4 


48 


/ 


t- 


f 


- ■ 


vn 


5 XX<9) =1000 
XX(13>=XKf t» 


SUBROUTINE  HAXHtRHO, FN0,P«,P£N,R1,Y,R2,T,I> 


DINEN?ION  X(?»  ,2r>) 


X(2)=8.o|E-i4 


T(2)=3.14£*t2 


101  c0RMAT (lx, *1  MAG IN ARY*) 


C2=FN0*YfI> 


^4*F£N*2tT) 


B=ri*C4-C2»33 


0=B*9-4.*A*C 


ft  1 I?  iw>Tir.w 


TPtO.UT.tr.O)  GO  TO  99 


W 1 


CHECKS  r 0?  »R3PE©  ROOTS 


T*fC3-C2*Y)  t (C2"Ci*Y) 


PPINT  10),Y,T,G 


R2«(C2-Y»C1I  /( t-Y) 


Y»Y*Y^l,i5E»y4 


G*G*5n.l&E»Q4 


99  PRINT  101 


